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1  INTRODUCTION 


This  projects  involves  the  development  of  methods  for  "High  Fidelity  Electronic 
Display  of  Digital  Mammographs” .  Cathodolumescent  medical  image  displays 
are  considered  including  specialized  cathode  ray  tubes  and  flat  panel  displays 
using  cold  cathode  emitters  and  microvacuum  cells  (fig.  1).  For  field  emitter 
display  devices,  a  small  cathode  generates  a  regulated  pulse  of  electrons  which 
is  accelerated  across  a  vacuum  gap  to  a  photoemissive  layer.  We  are  focusing 
in  this  project  on  the  visual  performance  of  the  photoemissive  layer. 

The  objective  of  our  research  is  to  establish  that  advanced  designs  for  the 
emissive  structure  of  cathodoluminescent  display  devices  can  achieve  the  perfor¬ 
mance  required  for  viewing  computer  images  obtained  from  digital  mammogra¬ 
phy  systems.  The  approach  being  investigated  is  to  building  the  photoemissive 
Conventional  CRT  devices  are  not  capable  of  meeting  the  stringent  display  re¬ 
quirements  of  medical  mammography  because  of  the  thick  glass  in  the  emissive 
structure.  We  think  that  a  microchannel  photoemissive  assembly  can  achieve 
the  needed  display  performance.  We  are  now  optimizing  the  design  of  the  pho¬ 
toemissive  assembly.  We  will  next  experimentally  establishing  the  ability  to 
achieve  high  brightness,  high  resolution  and  low  glare  with  a  wide  dynamic 
range. 


Figure  1:  Photoemissive  emissive  structure  model. 

In  most  display  devices,  image  degradation  can  be  traced  to  processes  asso¬ 
ciated  with  the  transport  of  light  from  the  point  of  generation  and  to  the  optical 
properties  of  the  emission  surface.  Transilluminated  film,  one  of  the  highest  fi¬ 
delity  display  devices  available,  modulates  the  light  intensity  by  diffusion  in  the 
layer  of  dense  silver  grains  in  the  emulsion.  This  layer  is  at  the  outer  surface 
and  separated  only  by  a  thin  overcoat  from  the  surface.  Consequently,  there  is 
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little  lateral  transport  of  light  which  would  degrade  contrast  or  resolution.  Ad¬ 
ditionally,  the  overcoat  has  a  rough  surface  causing  diffuse  scatter  of  ambient 
light  and  minimal  surface  glare. 

We  seek  a  photoemissive  assembly  which  generates  light  very  near  the  surface 
and  with  a  rough  outer  surface  so  as  to  appreciate  the  same  high  fidelity  as 
obtained  with  transilluminated  film.  The  general  approach  we  are  using  involves 
the  use  of  a  glass  microchannel  plate  with  an  overall  thickness  of  2  to  3  mm.  As 
illustrated  in  fig.  1,  electrons  from  the  cathode  elements  are  directed  into  the 
microchannels  to  a  thin  phosphor  layer  deposited  on  glass  with  a  thickness  of 
about  100  microns,  i.e.  similar  to  that  of  the  cover  plates  used  in  conventional 
microscope  slides.  Structural  rigidity  is  provided  by  the  microchannel  plate. 
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2  Body:  Year  One  Progress 

This  project  involves  four  technical  objectives: 

•  Model  the  photoemissive  performance  of  microchannel  assemblies. 

•  Fabricate  a  test  system. 

•  Design  and  construct  a  microchannel  emissive  structure. 

•  Perform  experimental  tests  and  compare  with  model  calculations. 

The  first  two  objectives  have  been  pursued  during  the  first  year  of  the  project  as 
was  planned  in  the  application.  Our  progress  in  these  two  areas  is  summarized 
in  the  following. 

2.1  Photoemissive  performance  modeling 

The  first  technical  objective  was  the  performance  modeling  of  the  proposed 
emissive  structure  assembly.  The  simulation  of  device  performance  relates  to 

3  different  aspects  which  differ  both  in  the  physical  processes  as  well  as  in  the 
applied  methods: 

•  Predicting  electron  trajectories  in  the  vacuum  cell, 

•  Predicting  electron  energy  deposition  in  the  emissive  layers,  and 

•  Simulation  of  the  light  transport  processes  in  the  emissive  structure. 

For  each  aspect  of  performance  modeling  we  describe  in  the  following,  the 
advancements  and  developments  and  results  made  during  this  first  half  of  the 
project  that  contributed  to  obtaining  significant  insight  on  the  performance  of 
the  proposed  emissive  structure  design. 

2.1.1  Modeling  electron  trajectories  in  vacuum  cells 


Methods 

To  study  the  effect  of  the  capillary  array  on  the  electron  trajectories  in  the 
proposed  design,  we  have  developed  a  simulation  code  (ELECTRA)  that  com¬ 
putes  the  solution  to  the  Poisson’s  equation  using  a  finite-difference  approach 
and  an  over-relaxation  technique.  For  the  analysis  of  the  impact  of  the  capillary 
structure,  the  full  3D  solution  may  be  obtained  using  a  cylindrical  coordinate 
system,  by  aligning  the  axis  of  symmetry  with  the  axis  of  the  capillary  tube, 
when  appropriate  boundary  conditions  are  used. 

An  important  aspect  of  computing  potential  solutions  in  vacuum  cells  for  the 
simulation  of  electron  trajectories  is  the  use  of  an  appropriate  mesh  size.  Large 
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array  sizes  provide  accuracy  while  increasing  the  computing  time  in  a  non-linear 
fashion  with  respect  to  array  dimensions.  To  achieve  good  solutions  for  the  elec¬ 
tric  potential,  we  implemented  a  3-stage  computation  scheme  (low-resolution, 
medium-resolution,  and  high-resolution).  By  computing  the  solution  at  each 
step,  and  using  a  bicubic  spline  technique  to  interpolate  the  expanded  solution, 
we  were  able  to  significantly  reduced  the  computing  time  for  a  given  solution, 
and  obtain  convergence  much  faster. 


Figure  2:  Isocontour  plot  for  a  region  close  to  the  focusing  plate. 


Figure  3:  Electron  trajectories  generated  with  ELECTRA  for  a  region  close  to 
the  focusing  plate. 

To  determine  when  the  iteration  has  converged  to  an  acceptable  solution, 
the  residual  at  a  given  point  in  the  solution  array  is  reported.  Typical  values  of 
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10-3  have  shown  to  be  acceptable  criteria  for  the  studied  cases. 

The  program  generates  electron  trajectories  using  the  Burlish-Stoer  method 
to  solve  the  equation  of  motion  with  adaptable  time  steps. [9]  The  results  are 
presented  in  isocontour  plots  and  trajectory  data  files  (fig.  2  and  3). 

Results 

We  have  applied  ELECTRA  to  study  the  effect  of  the  focusing  plate  electrode 
for  the  proposed  design.  For  a  3  mm  thick  structure,  a  1  mm  thick  microcapillary 
array  with  a  1  kV  total  acceleration  voltage  and  focusing  metal  electrode  at  V) 
was  considered.  The  metal  electrode  enters  the  channel  up  to  a  diameter  length, 
i.e.  50  /zm.  Electrons  are  assumed  to  be  emitted  from  the  cathode  forming 
a  parallel  extended  beam  for  normal  incidence  into  the  microcapillary  array. 
This  can  be  achieved  with  an  emitting  surface  such  as  a  field  emitter  cathode 
structure,  or  with  a  point  emitter  and  appropriate  electron  optics  components. 
The  goal  is  to  adjust  V)  in  order  to  alter  the  trajectories  to  direct  the  electrons 
into  the  tubes,  while  minimizing  the  impacts  on  the  focusing  electrode  plate 
and  on  lateral  tube  walls  (see  next  section  for  a  discussion  on  secondary  electron 
emission). 

Problems  encountered 

A  potential  problem  with  the  approach  being  used  involves  the  charging  of 
insulators  in  vacuum  under  electron  bombardment  which  has  been  studied  by 
several  authors.  [7,  6,  1]  When  an  energetic  electron  hits  an  insulator,  multiple 
secondary  electrons  are  generated.  A  fraction  of  them  are  directed  back  into 
the  vacuum.  This  fraction  is  called  ’’true”  secondary  electron  emission.  En¬ 
ergetic  backscattered  electrons  can  travel  significant  distances  away  from  the 
initial  landing  site.  The  low-energy  electrons  however,  can  in  some  scenarios, 
be  directed  back  against  the  surface,  creating  charge  accumulation  in  the  super¬ 
ficial  layers  of  the  insulator  material.  This  surface  charge  density  may  in  turn 
affect  the  local  electric  fields,  up  to  a  point  where  a  self-sustaining  equilibrium 
is  achieved.  In  this  case,  a  secondary  electron  emission  of  1  is  required.  When 
electric  fields  are  parallel  to  the  insulator  surface,  a  net  current  can  be  observed 
which  is  explained  by  a  hopping  transport  mechanism. [8] 

Inside  capillaries,  the  electric  field  is  normally  oriented  parallel  to  the  insu¬ 
lator  surface.  Two  different  cases  can  be  studied.  In  the  first,  an  energetic  elec¬ 
tron  may  hit  directly  the  interior  wall  of  the  tube,  before  reaching  the  phosphor 
layer.  Such  an  event  will  generate  multiple  less  energetic  secondary  electrons 
that  will  eventually  reach  the  anode.  This  condition  is  not  favorable  to  an  ef¬ 
ficient  electron  energy  to  luminescence  conversion,  and  thus  its  occurrence  has 
to  be  minimized  by  using  an  appropriate  voltage  Vf .  The  second  case  relates  to 
the  backscattered  electrons  coming  from  the  phosphor  material  that  may  hit  the 
surfaces  of  the  tube.  In  this  case,  the  wide  range  of  electron  energies  and  angles 
of  interaction  will  determine  the  ability  of  the  system  to  reach  the  equilibrium 
point. 
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Similar  scenarios  are  encountered  in  other  technologies  related  to  emissive 
displays.  For  instance,  Micron  Display  Technologies,  Inc.  has  recently  reported 
the  use  of  high  aspect  ratio  glass  posts  for  spacing  the  cathode  structure  from 
the  anode  layers.  [4]  The  use  of  a  thin  resistive  layer  coated  onto  the  inner  walls 
of  the  tubes  may  constitute  a  solution  to  this  problem,  although  a  mechanism  for 
controlling  the  resistivity  of  a  thin  layer,  as  well  as  a  suitable  coating  technique 
have  not  been  disclosed. 

This  problem  was  not  recognized  at  the  beginning  of  this  project.  At  the 
present  time,  we  are  addressing  the  issue  of  surface  charging  by  using  simulation 
methods  to  describe  the  energy  and  angular  distribution  of  secondary  electrons. 
We  will  also  determine  experimentally  under  what  circumstances  an  energetic 
beam  will  travel  through  the  tube  without  losses. 

2.1.2  Electron  energy  deposition  modeling 


Glare  from  backscattered  electrons 

Our  laboratory  previously  developed  an  electron/photon  transport  code, 
SKEPTIC  [5,  10],  specifically  intended  to  model  low  energy  electron  transport. 
We  have  used  SKEPTIC  to  characterize  the  electron  backscattering  from  the 
phosphor  layer  of  an  emissive  structure.  In  CRT  devices,  the  backscattered 
component  that  returns  into  the  vacuum,  may  hit  the  phosphor  layer  at  signif¬ 
icant  distances  from  the  original  landing  site,  and  generate  a  diffuse  signal.  We 
have  completed  computations  to  document  the  effect  of  this  backscattering  on 
image  glare.  For  this,  we  used  actual  layer  material  specifications  including  the 
aluminum  conductive  film  on  top  of  the  phosphor  grains,  and  the  conductive 
coatings  on  the  inside  of  the  tube.  We  have  used  a  typical  broad  band  ZnO:Zn 
phosphor  for  the  simulations.  Preliminary  results  show  that  the  tails  of  the  line- 
spread  function  for  incidence  at  the  center  will  not  affect  significantly  the  image 
quality  of  the  displayed  data  (fig.  4).  However,  in  color  CRTs,  grille  apertures 
or  shadow  masks  may  interfere  the  backscattered  electrons  and  generate  a  more 
important  and  local  degradation  effect. 

Electrons  backscattered  in  channels 

We  are  now  using  SKEPTIC  to  model  the  angular  and  energy  distribution  of 
backscattered  electrons  coming  from  the  interaction  of  energetic  electrons  with 
the  phosphor  in  the  proposed  emissive  structure  design.  This  will  contribute 
to  the  description  of  the  physics  of  electrons  traveling  inside  the  microcapillary 
array,  and  the  effect  of  backscattering  on  surface  charging. 

2.2  Modeling  of  light  transport  in  emissive  structures 

We  have  developed  a  Monte  Carlo  simulation  code  (DETECT-II)  to  study  the 
luminance  spread  functions  of  emissive  structures.  The  unique  features  of  the 
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Figure  4:  Line-spread  function  of  backscattered  electron  diffuse  signal  for  a 
typical  monochrome  CRT. 


code,  as  well  as  the  results  for  typical  CRT  and  flat  emissive  display  structures 
have  been  reported.  [2,  3]. 

In  addition,  simulation  of  degradation  of  image  quality  by  glare  using  2D 
Fast  Fourier  Transform  over  large  image  arrays  (2048  x  2048)  was  performed. 
The  application  of  this  method  to  digital  radiographs  was  described  recently. 
[2]  We  refer  to  the  proceeding  papers  which  are  attached  as  appendices  for  a 
complete  description  of  these  methods  and  results. 


Figure  5:  Test  device  assembly  for  testing  photoemissive  emissive  structure. 


2.2.1  Model  validation 

We  have  reported  the  glare  for  both  thick  and  thin  emissive  structures  in  the 
papers  noted  above.  We  are  now  performing  experiments  on  sections  of  thick 
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structures  extracted  from  medical  imaging  CRT  devices.  These  experiments 
will  establish  the  accuracy  of  our  models  as  well  as  document  the  performance 
deficit  of  current  CRT  devices. 

2.3  Fabrication  of  a  test  system 

2.3.1  Vacuum  subsystem 

The  test  system  for  emissive  structures  being  built  consists  of  a  vacuum  chamber 
with  a  cathode  at  one  end  and  an  emissive  structure  installed  at  the  other  end. 
Tests  are  planned  in  which  a  pattern  will  be  placed  on  the  back  surface  of  the 
emissive  structure  to  enable  resolution  and  glare  to  be  measured.  Used  in  the 
manner,  the  cathode  needs  only  to  produce  a  uniform  beam  of  electrons  with 
no  dynamic  pattern. 


detector 


Figure  6:  Cone-shaped  probe  for  glare  measurements. 

The  design  of  the  vacuum  system  has  been  completed  all  components  are  now 
being  installed  in  our  laboratory.  This  includes  a  low  vacuum  mechanical  pump, 
a  high  vacuum  turbopump,  high  and  low  vacuum  gauges,  and  an  experimental 
chamber.  The  configuration  of  the  vacuum  subsystem  is  detailed  in  an  attached 
appendix. 

2.3.2  Glare  performance  measurements 

Measurements  of  glare,  particularly  for  systems  with  low  glare,  require  measure¬ 
ment  systems  which  do  not  contribute  to  the  observed  glare.  The  measurements 
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we  proposed  were  to  all  be  done  with  a  CCD  observation  system.  This  will  still 
be  used  for  measurements  of  resolution  and  noise.  For  glare  measurements  we 
have  designed  an  optical  probe  to  be  used  with  a  laboratory  photometer.  The 
cone  shaped  luminance  probe  will  be  used  to  measure  luminance  spread  func¬ 
tion  tails  of  actual  emissive  structures  both  in  the  test  sytem  and  conventional 
display  devices  (fig.  6). 
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3 


CONCLUSIONS 


During  the  first  half  of  the  project  we  have  advanced  the  simulation  tools  to 
characterize  electron  and  optical  transport  in  emissive  structures.  By  using  the 
simulation  tools,  we  have  advanced  our  understanding  of  the  degrading  effects 
of  the  luminance  spread  in  thick  emissive  structures. 

In  the  next  half,  we  will  demonstrate  our  findings  on  improved  display  per¬ 
formance  of  such  an  emissive  structure  by  measuring  the  characteristics  of  a  test 
device  (fig.  5)  that  will  consist  of  a  parallel  extended  electron  beam  and  the 
proposed  emissive  structure  with  the  microcapillary  array  backing  the  anode 
phosphor  layer. 
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X:  Experimental  Chamber 

gl:  Active  Inverted  Magnetron  Gauge 
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Image  Degradation  by  Glare  in  Radiologic  Display  Devices 
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ABSTRACT 


No  electronic  devices  are  currently  available  that  can  display  digital  radiographs  without  loss  of  visual  informa¬ 
tion  compared  to  traditional  transilluminated  film.  Light  scattering  within  the  glass  faceplate  of  cathode-ray  tube 
(CRT)  devices  causes  excessive  glare  that  reduces  image  contrast.  This  glare,  along  with  ambient  light  reflection, 
has  been  recognized  as  a  significant  limitation  for  radiologic  applications.  Efforts  to  control  the  effect  of  glare  and 
ambient  light  reflection  in  CRTs  include  the  use  of  absorptive  glass  and  thin  film  coatings.  In  the  near  future, 
flat  panel  displays  (FPD)  with  thin  emissive  structures  should  provide  very  low  glare,  high  performance  devices. 
We  have  used  an  optical  Monte  Carlo  simulation  to  evaluate  the  effect  of  glare  on  image  quality  for  typical  CRT 
and  flat  panel  display  devices.  The  trade-off  between  display  brightness  and  image  contrast  is  described.  For 
CRT  systems,  achieving  good  glare  ratio  requires  a  reduction  of  brightness  to  30-40  %  of  the  maximum  potential 
brightness.  For  FPD  systems,  similar  glare  performance  can  be  achieved  while  maintaining  80  %  of  the  maximum 
potential  brightness. 

Keyword  List:  flat  panel  display,  emissive  display,  digital  radiography,  CRT,  glare. 


1  INTRODUCTION 


The  radiologic  information  acquired  by  current  digital  radiographic  systems  is  often  interpreted  using  film 
from  a  laser  printer.  The  electronic  display  of  digital  radiographs  will  require  devices  with  the  performance  of 
transilluminated  film.1  High  luminance,  low  noise  and  wide  dynamic  range  are  required,  as  summarized  in  Table 
I.  Computer  workstations  designed  for  displaying  and  interpreting  digital  radiographs  all  use  cathode-ray  tube 
monitors  (CRT).  However,  the  limitations  of  current  CRT  display  devices  as  compared  to  transilluminated  film 
have  been  recognized  and  include  issues  such  as  brightness,  resolution,  dynamic  range,  uniformity,  and  noise.2-6 

Typical  CRTs  require  a  thick  glass  panel  between  the  light  emission  sites  and  the  viewer.  Depending  on  the 
curvature  ratio  of  the  tube,  the  faceplate  thickness  may  range  between  1.3  and  2.5  cm.7  Multipath  light  scattering 
in  the  faceplate  adds  a  diffuse  background  (glare)  to  the  primary  signal  that  reduces  contrast.  The  degradation 


Total  dimensions 

35  cm  x  43  cm 

Pixel  dimensions 

160  n m  max.  /  80  /m l  nominal 

Array  size 

2000  x  2500  min.  /  4000  x  5000  nominal 

Peak  luminance 

2000  cd/m2 

Minimum  luminance 

5  cd/m2 

Color 

white 

Emission  distribution 

lambertian 

Noise  power  spectrum 

white 

S/N  for  each  pixel 

>  100 

S/N  for  .5  mm  region 

>  400 

Intrascene  dynamic  range 

400 

Greyscale 

1024  (log  or  perceptually  linear) 

Refresh  rate 

static  or  70  Hz 

Large  area  distortion 

1% 

Color  inversion 

none 

Table  1:  High  fidelity  display  requirements 


in  image  quality  is  significantly  more  severe  for  subtle  lesions  in  dark  regions  with  bright  surroundings.  To  reduce 
the  contrast  reduction  by  glare,  dark  tinted  glass  is  used  for  the  faceplate.  Transmission  can  be  as  low  as  13  % 
for  a  thickness  of  1.78  cm  (equivalent  to  an  absorption  coefficient  of  1.15  cm-1).8  On  the  other  hand,  flat  panel 
field-emission  displays9"13  have  the  potential  for  a  thin  faceplate  due  to  a  large  number  of  spacers  between  the 
electron  emitting  layer  and  the  phosphor.  In  this  paper,  the  low  frequency  degradation  of  image  quality  by  glare 
from  light  scattering  in  the  faceplate  of  emissive  display  devices  is  addressed.  The  effect  of  glare  for  both  CRT  and 
flat  panel  display  emissive  structures  is  computationally  modeled  for  test  pattern  images.  Finally,  the  influence 
of  glare-reducing  glass  absorbers  on  both  glare  and  display  brightness  is  reported. 


2  GLARE  IN  EMISSIVE  DISPLAYS 


Computational  models  have  been  used  to  simulate  the  effect  of  device  characteristics  on  display  performance 
including  the  frequency  response  of  electronics,  electron  beam  spot  size,  and  phosphor  granularity.14,15  Com¬ 
putational  modeling  can  also  be  used  to  simulate  glare.  This  requires  accurate  knowledge  of  all  processes  that 
propagate  light  from  bright  spots  in  an  image  to  points  at  substantial  distance,  i.e.  the  tails  of  the  device  point 
spread  function  (PSF).  This  work  has  been  motivated  by  a  desire  to  obtain  accurates  PSFs  to  model  glare.  While 
the  tails  of  these  PSFs  can  have  very  low  values,  they  can  contribute  significant  signal  because  of  their  spatial 
extent. 

Glare  measurements  on  actual  CRT  emissive  structures  have  been  recently  reported.16  These  glare  measure¬ 
ments  and  glare  measurement  results  are  strongly  dependent  on  the  nature  of  the  test  image.  However,  no 
standard  metric  has  been  adopted  for  reporting  the  glare  characteristics  of  actual  display  devices.  In  this  study, 
a  test  pattern  that  we  previously  described  was  used.17  The  pattern  consists  of  an  inner  dark  circular  spot,  sur¬ 
rounded  by  a  bright  outer  circle  with  a  ratio  of  radii  of  10.  The  rest  of  the  scene  is  kept  at  the  same  dark  level  of 
the  inner  spot  (see  Fig.  1).  The  glare  ratio  is  defined  as  the  difference  in  display  luminance  between  the  central 
pixel  in  the  pattern  and  the  brightest  pixel  in  the  bright  region. 

The  impact  of  electron  backscattering  on  display  performance  has  not  been  considered  in  previous  studies. 
For  an  A1  layer,  the  backscattering  fraction  for  an  electron  beam  energy  of  10  keV  is  about  12.5-15  %  and 
14-15  %  for  30  keV. 18,19  The  backscattered  electrons  are  subjected  to  the  electromagnetic  fields  present  in  the 


vacuum  region.  The  backscattered  electrons  will  eventually  impinge  on  the  phosphor  at  a  distance  from  its  initial 
interaction  position  which  is  a  function  of  the  backscattering  angle  and  energy.  The  backscattering  energy,  which 
is  a  function  of  the  backscattering  angle,  will  eventually  result  in  the  backscattered  electron  impinging  into  the 
phosphor  at  a  certain  distance  from  its  initial  interaction  position.  This  results  in  an  electronic  contribution  to 
glare  that  is  not  considered  in  this  work. 


Figure  1:  Test  pattern  for  glare  measurements. 


3  MODELING  GLARE  IN  DISPLAY  DEVICES 


3.1  Monte  Carlo  calculation  of  the  point  spread  function 

In  a  cathodoluminescent  display  device,  the  emissive  structure  includes  all  components  that  convert  energetic 
electrons  to  visible  light.  Typical  emissive  structures  components  include:  cathodoluminescent  phosphors,  con¬ 
ductive  layers,  reflective  or  absorptive  films,  transparent  support  (faceplate),  and  anti-reflective  and  anti-glare 
materials.  In  this  work,  we  have  modeled  the  glare  characteristics  of  2  types  of  emissive  structures.  All  surfaces 
are  modeled  as  perfectly  flat.  The  angular  distribution  emerging  from  the  phosphor  layer  was  assumed  to  be 
lambertian. 

Structures  CRT-A  and  CRT-B  are  typical  CRT  display  devices.  Both  have  a  faceplate  thickness  of  1.6  cm. 
CRT-A  has  an  A1  conductive  and  reflective  layer  on  the  back  of  the  phosphor  layer  producing  a  reflection  of  90 
%.  CRT-B  has  a  perfectly  absorbing  black  matrix  coated  in  between  the  phosphor  dots,  with  an  aperture  ratio 
of  50  %.  For  each  CRT  structure,  5  levels  of  absorption  coefficient  were  considered  (0.0,  0.2,  0.4,  0.6,  and  1.0 
cm-1).  Emissive  structures  FPD-A  and  FPD-B  are  typical  of  flat-panel  emissive  displays.  The  thicknesses  are 
0.3  and  0.1  cm  respectively.  For  each  FPD  structure,  8  levels  of  absorption  coefficient  were  considered  (0.0,  0.2, 
0.4,  0.6,  1.0,  3.0,  5.0  and  9.0  cm-1).  Cross-sections  of  the  emissive  structures  are  presented  in  figures  2  and  3. 

In  previous  work,  we  described  Detect-II,  an  optical  Monte  Carlo  code,  capable  of  simulating  the  optical 
transport  processes  in  emissive  structures.  Unique  features  of  the  code  include  special  binning  for  simulating 
display  performance  characteristics,  and  photon  polarization  tracking.17  Detect-II  was  used  to  obtain  the  PSFs 
of  different  emissive  structures.  The  results  of  the  Monte  Carlo  calculation  yield  the  number  of  photons,  F0,  that 
originate  from  the  source  point  and  are  observed  within  the  solid  angle  of  the  observer.  Secondly,  the  number 


A1  backing  with  90  %  diffuse  reflectance  Black  matrix  with  50  %  aperture  ratio 


(a)  (b) 

Figure  2:  Cross-sections  of  emissive  structures  CRT-A  (a)  ,  and  CRT-B  (b). 


50  %  diffuse  reflectance  50  %  diffuse  reflectance 


(a)  (b) 

Figure  3:  Cross-sections  of  emissive  structures  FPD-A  (a),  and  FPD-B  (b). 


of  photons  observed  within  the  solid  angle  of  the  observer  but  last  associated  with  a  point  on  the  emission  layer 
which  is  located  at  a  radius  from  r  to  Ar  from  the  source  point  are  binned  within  N  bins  out  to  a  distance  NAr(P* 
for  i  =  1..1V).  This  is  then  converted  to  the  differential  probability  of  light  emission  per  unit  area, 


L(i)  = 


r  (  *\ 

L(,)  =  iA 

Pi/Po 


;r[((i  +  l)Ar)2  —  (iAr)2] 


for, 


r  =  iAr 
i  =  1..N 

where  L(i )  is  a  discrete  representation  of  the  continuous  luminance  point  spread  function  L(r).  L(r)  can  be 
interpreted  as  the  observed  luminance  in  cd/m2  resulting  from  a  point  source  located  at  the  origin  of  the  emissive 
surface  with  a  luminous  intensity  of  1  cd. 


Because  7(0)  is  associated  with  a  point  source  and  goes  to  a  very  large  number  when  small  Ar  are  used  for  the 
binning  of  the  Monte  Carlo  code,  we  define  L(i)  only  for  values  of  i  ranging  from  1  to  N.  Instead,  the  luminous 
intensity  of  the  central  peak,  7(0),  was  defined  as  the  fraction  of  unscattered  light  that  reachs  the  viewer  from 
a  point  source  producing  a  luminous  intensity  of  1  cd.  1(0)  describes  the  brightness  of  images  not  degraded  by 
glare.  For  a  digital  display,  the  luminance  associated  with  an  emissive  structure  having  an  apparent  luminous 
intensity  of  1(0)  per  pixel  is  7(0)/pp,  where  pp  is  the  number  of  pixels  per  m2.  1(0)  can  be  deduced  from  Pq 
and  the  lumens  of  light  associated  with  the  number  of  histories  started  at  the  phosphor  layer.  For  the  emissive 
structures  studied  in  this  work,  Pq  is  the  same  in  the  absence  of  absorbers.  We  therefore  have  normalized  1(0)  to 
1000.  The  relative  reduction  of  7(0)  with  absorption  is  used  as  a  measure  of  brightness. 
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Figure  6:  Processing  steps  for  convolution  of  radiographic  images  with  glare  PSF. 


Monte  Carlo  data  was  computed  out  to  a  distance  of  20  cm  for  100  bins  with  a  spacing  of  2  mm.  A  two 
dimensional  point  spread  function  ( L(x ,  y))  was  interpolated  from  this  for  a  20  x  20  cm  area  with  2048  x  2048 
pixels  having  dimensions  of  0.2  x  0.2  mm.  All  points  out  to  a  radius  of  2  mm  were  set  equal  to  L(  1).  The 
remaining  p'oints  at  larger  radius  were  linearly  interpolated  from  L(i).  All  points  other  that  the  center  point  were 
then  converted  from  cd/m2  to  cd/pixel  by  multiplying  each  value  by  an  area  per  pixel  of  0.04  mm2  to  correspond 
with  digital  images  having  a  0.2  mm  pixel  size.  Finally,  the  point  spread  function  was  normalized  so  that  the 
sum  of  all  points  in  the  2-dimensional  discrete  array  was  1.0.  This  insures  that  the  convolution  with  the  PSF 
preserves  the  average  brightness  of  a  displayed  scene.  Implicitly,  the  display  of  an  image  with  high  glare  is  thus 
assumed  to  be  adjusted  to  preserve  the  same  average  brightness. 

Figures  4  and  5  show  the  Pi/  Pq  of  the  4  emissive  structures  computed  in  this  work.  The  amplitude  of  the 
tails  of  these  functions  decreases  with  the  absorption  level.  For  CRT-A,  the  tail  extends  up  to  a  distance  of  20 
cm,  with  a  large  amplitude  for  low  absorption  levels.  For  CRT-B,  the  overall  same  profile  is  observed,  although 
the  tail  amplitudes  are  lower  and  decrease  more  rapidly  with  distance  as  compared  to  the  PSFs  of  CRT-A.  For 
the  flat-panel  emissive  structures,  tails  have  a  greater  amplitude  close  to  the  origin  and  a  rapid  fall-off.  For  the 
thinner  structure  (FPD-B),  the  amplitudes  of  the  tails  are  still  significant  close  to  the  origin  but  rapidly  decay 
to  very  low  numbers  with  distance. 


3.2  Computational  simulation  of  glare  for  displayed  radiographs 

In  this  study,  we  assumed  that  the  display  of  radiographs  in  electronic  devices  can  be  modeled  by  convolving 
radiographic  images  in  luminance  units  with  the  2-dimensional  PSF,  L(x,  y ).  The  convolution  process  is  performed 
on  images  in  luminance  units,  since  the  degrading  effect  of  glare  is  linear  after  the  visible  light  has  been  generated  in 
the  phosphor  layer.  For  simulation  purposes,  images  in  film  density  units  are  first  obtained  by  (a)  computationally 
generating  a  test  pattern  image  of  film  density,  (b)  digitizing  available  radiographic  films,  or  (c)  converting  CR 


(a)  (b) 

Figure  7:  Test  pattern  image  center  data  row  for  different  absorption  coefficients  in  cm-1  after  convolution  with 
PSFs  of  CRT-A  (a),  and  CRT-B  (b). 
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Figure  8:  Test  pattern  image  center  data  row  for  different  absorption  coefficients  in  cm  1  after  convolution  with 
PSFs  of  FPD-A  (a),  and  FPD-B  (b). 
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Figure  9:  Glare  ratios  and  central  pixel  relative  brightness  as  a  function  of  absorption  coefficient  for  (a)  CRT, 
and  (b)  flat -panel  emissive  structures. 


data  in  units  of  log(signal)  to  film  density  units  using  parameterized  HD  curves.  These  images  were  then  converted 
from  density  to  luminance  units  using  the  following  relationship:  L  =  Lmax  *  10~FD ,  where  Lmax  is  the  viewbox 
luminance  (assumed  at  3000  cd/m2),  and  FD  is  the  film  density. 

To  simulate  the  degradation  of  image  quality  by  glare,  2K  x  2K  images  were  convolved  in  the  frequency 
domain  with  2K  x  2K  filters  generated  from  L(x,  y).  Because  L(x,  y)  has  extended  tails,  L(x,  y)  arrays  of  4Iv  x 
4K  are  needed  for  exact  solutions.  Our  use  of  2K  x  2K  kernels  may  introduce  aliasing  artifacts,  although  none 
were  observed  in  the  images  studied  to  date. 

L(x,  y)  was  first  transformed  using  a  2-dimensional  fast  Fourier  transform  (FFT).  Since  this  function  is  sym¬ 
metric,  we  stored  the  magnitude  of  the  Fourier  coefficients  in  real  arrays  that  were  then  used  as  filter  functions. 
The  2K  luminance  images  were  read,  Fourier  transformed  and  filtered  by  multiplying  the  complex  transformed 
image  by  the  filter  value.  The  complex  array  resulting  from  this  multiplication  was  then  converted  back  to  lumi¬ 
nance  space  by  an  inverse  FFT  procedure.  After  the  convolution  process  was  performed,  images  were  converted 
back  to  film  density  units.  All  computations  were  performed  on  a  DEC  Alphasystem  2000.  The  convolution  pro¬ 
cess  for  the  2K  images  including  reading  and  writing  operations,  conversion  to  luminance  units,  FFT,  filtering, 
inverse  FFT  and  conversion  to  film  density  units  was  coded  in  fortran90  and  required  about  60  s  per  image.  The 
processing  steps  are  summarized  in  Fig.  6.  Provision  was  made  to  extract  a  subregion  in  film  density  units  form 
the  undegraded  and  glare  degraded  image,  i.e.  A  and  B  in  Fig.  6.  This  was  done  to  permit  observer  studies 
where  A  and  B  are  printed  on  film  using  a  laser  printer. 

Using  this  method,  circular  test  patterns  of  the  type  shown  in  Fig.  1  were  convolved  with  the  filters  computed 
for  the  4  emissive  structures.  Figures  7  and  8  show  the  results  of  computed  film  density  along  a  row  through 
the  center  of  a  test  pattern  with  a  10  cm  diameter  bright  circle  with  a  1  cm  diameter  dark  center.  A  diffuse 
glare  component  can  be  seen  surrounding  the  bright  circle  with  the  amount  of  glare  being  inversely  related  to 
the  amount  of  absorber.  The  glare  ratio  was  computed  as  the  ratio  of  the  luminance  in  the  bright  area  to  the 
luminance  in  the  dark  spot.  Figure  9  shows  the  reduction  in  glare  and  improvement  in  glare  ratio  that  can 
be  achieved  by  increasing  the  amount  of  absorbers.  Also  shown  is  the  decrease  in  overall  brightness  caused  by 
absorption. 


glare  ratio 

Figure  10:  Relative  brightness  as  a  function  of  glare  ratio  for  emissive  structures  CRT-A,  CRT-B,  FPD-A,  and 
FPD-B. 


4  DISCUSSION 

In  this  study  we  computed  results  for  circular  test  patterns  of  4  cm  diameter  and  20  cm  diameter  in  addition 
to  the  reported  results  for  10  cm  diameter.  The  glare  ratios  measured  for  the  different  patterns  are  different 
because  the  component  in  the  dark  center  is  integrated  over  a  larger  radius.  Further  work  should  be  done  to 
identify  which  test  pattern  is  best  with  respect  to  predicting  adverse  glare  degradation  in  clinical  radiographs. 

We  have  conducted  preliminary  studies  using  clinical  radiographs  obtained  from  a  digital  radiography  system. 
For  both  chest  and  skeletal  radiographs,  small  lesions  were  added  to  the  log(signal)  data  in  regions  where  glare 
degradation  was  anticipated.  Subjective  observation  of  this  data  indicates  that  a  glare  ratio  of  250  measured  with 
a  10  cm  diameter'test  pattern  does  not  create  appreciable  degradation  in  the  clinical  image.  Further  work  will 
be  necessary  to  understand  how  differing  amounts  of  glare  impact  radiologic  observation. 

From  our  results  it  is  apparent  that  emissive  structures  which  are  designed  to  minimize  glare  will  have  reduced 
brightness.  This  is  illustrated  in  Fig.  10  where  we  have  plotted  the  relative  brightness,  1(0),  as  a  function  of 
the  glare  ratio  for  all  cases  considered.  For  the  same  glare  ratio,  the  black  matrix  structure  is  seen  to  perform 
better  than  a  conventional  monochrome  CRT  although  in  this  simulation  we  have  not  accounted  in  differences 
of  phosphor  luminous  efficiency.  Notably,  a  significant  improvement  is  seen  in  very  thin  flat  panel  displays  for 
which  the  high  frequency  of  multipath  reflections  allow  the  phosphor  to  damp  the  tail  of  the  PSF. 
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Abstract 

The  luminance  spread  function  of  thin  emissive  dis¬ 
plays  has  been  modeled  using  an  optical  Monte  Carlo 
code  to  describe  light  transport.  Line-spread  and  mod¬ 
ulation  transfer  functions  have  been  computed  for  differ¬ 
ent  emissive  structure  thicknesses."  Simulation  of  image 
degradation  performed  by  a  convolution  process  showed 
that  even  for  very  thin  faceplates,  control  of  lateral  light 
diffusion  is  needed  to  minimize  the  contrast  degradation 
by  glare.  When  compared  to  thick  emissive  structures 
typical  of  cathode-ray  tubes,  thin  structures  are  found 
to  be  capable  of  high  quality  due  to  more  frequent  ab¬ 
sorption  by  the  phosphor  layer. 


For  some  applications  such  as  graphic  arts  and  medi¬ 
cal  imaging,  display  devices  with  very  good  image  quality 
are  required.  High  brightness,  low  noise,  wide  luminance 
range  and  low  glare  are  specifically  required  for  digital  ra¬ 
diographic  display  devices.5  To  depict  good  contrast  over 
the  wide  range  of  signals  in  a  digital  radiograph  requires 
a  dynamic  range  of  100.  To  minimize  the  degradation  of 
contrast  in  regions  of  minimum  luminance,  the  glare  from 
surrounding  bright  regions  should  be  less  than  0.25  of  the 
minimum  luminance  (i.e.,  1/400  of  the  maximum  lumi¬ 
nance).  This  work  examines  the  design  requirements  for 
achieving  this  low  glare  in  thin  emissive  display  devices. 

Currently,  workstations  designed  for  displaying  and  in¬ 
terpreting  digital  radiographs  all  use  cathode  ray  tube 
devices  (CRTs).  Multipath  light  scattering  in  CRT  face¬ 
plates  causes  extended  tails  in  the  luminance  spread  func¬ 
tion.6,1  Black  matrix  coatings  and  glass  absorption  are 
approaches  that  have  been  used  to  control  glare  from 
optical  scattering  at  the  expense  of  display  brightness. 
In  addition  to  the  optical  transport  processes,  electron 
backscattering  in  the  vacuum  tube  has  been  reported  to 
contribute  to  glare.4,10  In  general,  the  image  quality  of 
CRT  devices  has  not  been  adequate  for  displaying  digital 
radiographs  in  medical  diagnostic  applications. 

On  the  other  hand,  flat  emissive  displays  have  the 
potential  for  better  image  quality  than  is  possible  with 
a  CRT.  For  these  devices  a  thin  faceplate  can  be  em¬ 
ployed  due  to  the  large  number  of  spacers  between  the 


electron  emitting  layer  and  the  phosphor.  Extended  elec¬ 
tron  backscattering  is  not  present  and  lateral  transport 
of  optical  photons  can  be  controlled  more  effectively.  For 
thin  emissive  structures,  no  results  have  been  reported  on 
the  optical  transfer  characteristics  and  its  effect  on  glare, 
brightness,  and  resolution.  In  this  paper,  we  examine 
the  optical  transport  processes  in  the  faceplate  of  thin 
emissive  display  devices  and  its  effect  on  low  frequency 
image  signals.  The  effect  of  glare  for  flat  panel  emissive 
structures  is  computationally  modeled  and  the  character¬ 
istics  of  the  luminance  spread  functions  are  reported.  In 
addition,  the  influence  of  glass  absorption  on  glare  and 
display  brightness  is  analyzed. 

The  point  spread  function  and  modulation  transfer 
function  (MTF)  are  commonly  used  to  describe  the  res¬ 
olution  of  imaging  devices.  The  glare  in  a  display  device 
may  also  be  described  by  the  point  spread  function  if 
care  is  taken  to  document  the  low  signals  in  the  tail  of 
the  function  which  can  extend  for  considerable  distances. 
These  tails  are  associated  with  a  low  frequency  drop  in 
the  MTF  which  can  be  used  to  describe  contrast  reduc¬ 
tion  due  to  glare.  We  have  used  numeric  Monte  Carlo 
methods  in  this  work  to  deduce  the  2-dimensional  (2D) 
point  spread  function  of  specific  systems  over  distances 
equal  to  the  full  display  size.  This  is  then  used  to  evalu¬ 
ate  measures  of  image  quality. 

II.  Computational  methods 
A.  Luminance  spread,  function 

The  modeling  of  the  light  transport  processes  that  oc¬ 
cur  in  emissive  structures  can  be  done  using  either  ray 
tracing  or  Monte  Carlo  methods.  Ray  tracing  techniques 
rely  upon  a  priori  knowledge  of  the  system  response  to  de¬ 
termine  the  principal  ray  paths  that  will  contribute  to  the 
output.  On  the  other  hand,  Monte  Carlo  methods  track 
the  trajectories  of  a  large  number  of  optical  photons  to  de¬ 
scribe  the  statistically  averaged  output.  Absorption  and 
scattering  processes  are  modeled  with  the  Monte  Carlo 
method  by  randomly  sampling  the  distance  to  each  event 
from  probability  distribution  functions.  Absorption  and 
scattering  are  described  by  linear  interaction  coefficients 
(cm-1).  Monte  Carlo  methods  can  handle  a  wide  variety 
of  physical  cases  and  are  computationally  practical  when 
using  high  speed  computers. 


I.  Introduction 
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To  model  the  luminance  spread  functions  in  emissive 
displays,  we  have  developed  a  Monte  Carlo  optical  trans¬ 
port  code  written  in  fortran  90  (Detect-II).  The  geome¬ 
try  is  described  in  3D  cartesian  coordinates  using  slabs 
with  orthogonal  division  into  volumetric  cells.  The  sur¬ 
face  definitions  include  rough  surfaces,  mirrors,  perfect 
absorbers,  partial  diffuse  absorbers  with  lambertian  emis¬ 
sion  and  thin  films.  Photon  histories  are  started  as  point 
or  planar  sources  with  lambertian  or  isotropic  angular 
distribution.  Individual  photons  are  allowed  to  be  ab¬ 
sorbed  or  scattered  in  the  medium.  A  unique  feature  of 
Detect-II  is  the  tracking  of  photon  polarization.  In  very 
thin  structures,  multiple  light  reflections  occur,  and  po¬ 
larization  dependencies  in  the  Fresnel  coefficients  have  to 
be  considered. 

Neglecting  the  effect  of  finite  wavelength,  the  light 
transport  is  simulated  by  tracking  individual  photons  us¬ 
ing  a  physical  optics  description.  Reflection  and  trans¬ 
mission  coefficients  given  by  Fresnel’s  equations  are  then 
interpreted  as  probabilities.  The  transmission  and  reflec¬ 
tion  probabilities  for  parallel  and  perpendicular  compo¬ 
nents  are:3 

_  sin{26i)sin{26-z) 

sin2  {6 1  +  62)cos2{6\  —  62) 

sin(2^i)sm(2^) 

X,T  sin2  (6 1  +62) 

tan2(9i—92) 
tan2(6i  +  02) 

p  -  gz’n2(fli  ~  ftp 
1,'R  sin2(9i+92) 

where  9\  is  the  incidence  angle  and  62  is  given  by  Snell’s 
law. 

An  important  aspect  of  the  code  is  a  method  for  bin¬ 
ning  results  which  is  designed  for  display  modeling.  To 
simulate  display  performance,  the  viewer  is  assumed  to 
observe  the  image  from  a  direction  normal  to  the  surface. 
Photons  are  tracked  from  the  source  to  the  front  surface 
of  the  display  and  those  within  a  finite  solid  angle  about 
the  normal  are  binned  into  discrete  regions  on  the  surface. 
A  solid  angle  corresponding  to  a  6°  cone  is  used  rather 
than  the  typical  solid  angle  of  a  human  observer  in  order 
to  make  efficient  use  of  the  photon  histories.  The  bin¬ 
ning  of  the  emerging  position  is  done  by  backprojecting 
the  light  path  to  a  virtual  plane  of  emission.  This  ac¬ 
counts  for  the  difference  in  index  of  refractions  between 
air  and  glass.  This  allows  binning  from  a  large  solid  an¬ 
gle  to  a  common  virtual  focal  plane  inside  the  emissive 
structure. 


To  avoid  correlations  in  multidimensional  space  when 
a  large  number  of  histories  are  needed,  a  portable,  very- 
long  period  (2144)  lagged  Fibonacci  random  number  gen¬ 
erator  is  used.8,7 


50  %  diffuse  reflectance 


Fig.  1.  Cross-section  of  thin  emissive  structure  model. 

In  a  thin  cathodoluminescent  display  device,  the  emis¬ 
sive  structure  includes  all  components  that  convert  ener¬ 
getic  electrons  to  visible  light.  Typical  emissive  struc¬ 
tures  components  include:  cathodoluminescent  phos¬ 
phors,  conductive  layers,  reflective  or  absorptive  films, 
transparent  support  (faceplate),  and  antireflective  and 
antiglare  coatings.  In  this  work,  we  have  modeled  the 
glare  characteristics  of  thin  emissive  structures  with  dif¬ 
ferent  faceplate  thicknesses  (0.5,  1.0,  2.0,  and  3.0  mm) 
and  for  absorption  levels  (0.0,  1.0,  3.0,  5.0,  9.0,  and 
12.0  cm-1).  All  surfaces  are  modeled  as  perfectly  flat. 
Light  entering  the  phosphor  layer  is  assumed  to  return 
into  the  faceplate  with  a  probability  of  0.5  with  a  lamber¬ 
tian  angular  distribution.  A  cross-section  of  the  emissive 
structure  is  presented  in  figure  1. 

In  previous  work,  we  described  methods  to  obtain  the 
2D  luminance  spread  function  from  the  Monte  Carlo  re¬ 
sults.1  A  brief  summary  is  given  in  this  paragraph  to  fa¬ 
cilitate  the  interpretation  of  the  results  presented  in  this 
paper.  The  Monte  Carlo  calculation  yields  the  number 
of  photons,  Pq,  that  originate  from  the  source  point  and 
are  observed  within  the  solid  angle  of  the  observer  asso¬ 
ciated  with  a  point  on  the  emission  layer  which  is  located 
at  a  radius  from  r  to  Ar  from  the  source  point.  This  is 
then  converted  to  the  discrete  differential  probability  of 
light  emission  per  unit  area  L(i),  which  is  a  discrete  rep¬ 
resentation  of  the  continuous  luminance  spread  function 
L(r).  L(r )  can  be  interpreted  as  the  observed  luminance 
in  cd/m2  resulting  from  a  point  source  located  at  the  ori¬ 
gin  of  the  emissive  surface  with  a  luminous  intensity  of 
1  cd.  For  the  central  peak,  1(0),  is  defined  as  the  frac¬ 
tion  of  unscattered  light  that  reachs  the  viewer  from  a 
point  source  producing  a  luminous  intensity  of  1  cd.  1(0) 
describes  the  brightness  of  images  not  degraded  by  glare. 

B.  Image  processing 

We  use  the  2D  luminance  spread  function  to  simulate 
image  degradation  by  convolving  this  function  with  val¬ 
ues  for  an  image  scene.  These  convolutions  are  done  in 
the  frequency  domain  using  a  2D  Fourier  transform.  To 
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describe  the  degradation  both  the  line-spread  function 
(LSF)  and  the  MTF  are  deduced.  The  ID  luminance 
LSF  is  computed  by  projecting  the  2D  luminance  spread 
function  and  normalizing  the  ID  vector  to  unity.  A  ID 
Fourier  transform  is  then  performed  on  the  LSF  to  obtain 
the  MTF. 

Glare  measurements  on  CRT  emissive  structures  are 
commonly  performed  with  a  large  bright  region  and  a 
small  black  spot.9,4,10  The  ratio  of  luminances  in  the 
bright  field  and  the  dark  spot  is  defined  as  the  glare  ra¬ 
tio.  The  results  of  these  glare  measurements  are  strongly 
dependent  on  the  nature  of  the  test  image.  However,  no 
standard  metric  has  been  adopted  for  reporting  the  glare 
characteristics  of  actual  display  devices.  In  this  study, 
we  use  a  circular  test  pattern  that  we  have  previously  de¬ 
scribed.2  The  pattern  consists  of  an  inner  dark  circular 
spot,  surrounded  by  a  bright  outer  circle  with  a  ratio  of 
radii  of  10.  The  rest  of  the  scene  is  kept  at  the  same  dark 
level  of  the  inner  spot.  The  glare  ratio  is  defined  as  the 
difference  in  display  luminance  between  the  central  pixel 
in  the  pattern  and  the  brightest  pixel  in  the  bright  region. 
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Fig.  2.  Processing  steps  for  simulation  of  image  quality 
degradation  by  glare  by  convolution  with  luminance 
spread  functions 


A  high  fidelity  display  device  should  replicate  the  ap¬ 
pearance  of  transilluminated  film.  To  model  this,  we  as¬ 
sume  that  the  display  of  radiographs  in  an  electronic  de¬ 
vice  can  be  modeled  by  convolving  radiographic  images 
in  luminance  units  with  an  interpolated  2D  PSF.  Glare 
pattern  images  are  first  converted  from  film  density  to 
luminance  units  using  the  following  relationship: 

L  —  Lmax  *  10~FD 


where  Lmax  is  the  viewbox  luminance  (assumed  at 
3000  cd/m2),  and  FD  is  the  film  density.  The  convo¬ 
lution  process  is  then  performed  on  images  in  luminance 
units,  since  the  degrading  effect  of  glare  is  linear  after  the 
visible  light  has  been  generated  in  the  phosphor  layer  (see 
figure  2).  To  simulate  the  degradation  of  image  quality 
by  glare,  2 K  x  2 K  images  were  convolved  in  the  frequency 
domain  with  2 1\  x  2A'  filters  generated  from  the  2D  PSF.1 

III.  Results  and  Discussion 

Figure  3  shows  the  LSF  of  emissive  structures  consid¬ 
ered  in  this  work  with  no  glass  absorption.  For  thin  emis¬ 
sive  structures  with  faceplate  thickness  less  than  3  mm, 
the  LSFs  show  high  amplitudes  close  to  the  origin  with 
tails  that  decay  rapidly  to  very  low  numbers.  Moreover, 
the  LSF  for  0.5  and  1  mm  have  similar  features,  differing 
only  at  large  distances  and  very  close  to  the  origin.  At 
3  mm,  a  peak  is  observed  that  defines  a  halo  around  the 
center  spot.  The  amplitudes  at  short  distances  are  lower 
than  for  thinner  faceplates  by  a  factor  of  at  least  0.25.  In 
figure  4,  the  MTFs  of  the  same  emissive  structures  are 
plotted.  The  drop  in  amplitude  is  approximate  constant 
for  all  frequencies,  up  to  the  Nyquist  limit. 


Fig.  3.  Line-spread  functions  for  different  faceplate 
thicknesses  with  no  glass  absorption  (for  a  pixel  size 
of  0.2  mm). 


To  determine  if  an  emissive  structure  is  capable  of 
achieving  the  high  glare  ratio  needed  for  radiologic  inter¬ 
pretation  (i.e. ,  about  400),  we  have  computed  test  pattern 
glare  ratios  for  emissive  structures  of  different  thicknesses 
and  with  different  absorption  levels. 

Figure  5  illustrates  the  relative  brightness  1(0)  as  a 
function  of  the  glare  ratio  for  the  emissive  structures 
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Fig.  4.  Modulation  transfer  functions  for  different  face¬ 
plate  thicknesses  with  no  glass  absorption  (for  a  pixel 
size  of  0.2  mm). 


Fig.  5.  Relative  brightness  as  a  function  of  the  glare  ra¬ 
tio  for  different  faceplate  thicknesses.  For  each  curve, 
glass  absorption  is  varied  to  decrease  J(0)  and  in¬ 
crease  the  glare  ratio. 


considered.  Higher  glare  ratios  are  possible  to  achieve 
by  increasing  the  absorption  level,  although  decreasing 
brightness.  To  obtain  a  given  glare  ratio,  an  increase  in 
faceplate  thickness  results  in  a  reduction  of  the  display 
brightness.  For  a  glare  ratio  of  400,  the  brightness  has 
to  be  reduced  by  37.5  %  if  the  thickness  of  the  emissive 
structure  is  increased  from  1  to  3  mm.  Again,  similar 
characteristics  are  observed  for  faceplate  thicknesses  of 
0.5  and  1  mm.  As  seen  in  figure  3,  the  similar  features 
in  the  LSFs  may  be  responsible  for  the  similarities  in  the 
achievable  brightness  for  a  given  glare  ratio. 


IV.  Summary 

The  luminance  spread  functions  of  thin  emissive  struc¬ 
tures  typical  of  fiat  cathodoluminescent  displays  have 
been  modeled  with  an  optical  Monte  Carlo  code.  The 
tails  of  the  spread  functions  show  that  even  for  thin  struc¬ 
tures,  the  light  scattering  within  these  structures  has  to 
be  controlled  in  order  to  achieve  glare  ratios  suitable  for 
high  performance  applications.  The  use  of  absorption  in 
the  faceplate  has  been  shown  to  be  a  possible  solution,  al¬ 
though  affecting  the  display  brightness.  Other  solutions 
may  include  the  use  of  a  black  matrix  coating  in  between 
the  phosphor  dots. 
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